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Abstract: Commercial farming of anuran species that are declining in 
the wild raises a need to discriminate wild from farmed frogs. We hy¬ 
pothesized wild frogs might have extended hindlimbs due to greater 
frequency or intensity of jumping relative to farmed frogs, highlighting a 
morphometric approach to discrimination of wild from farmed frogs 
using hindlimb length. In the present study, Dybowski’s frog (Rana dy¬ 
bowskii) was used to test this hypothesis. We measured body mass (M b ) 
and hindlimb length (L h ) of 2-year old farmed frogs and wild frogs aged 
2 to 5 years. Dybowski’s frog demonstrated significant dimorphism in 
M b and L h . M b was significantly greater among farmed 2-year old frogs 
in both sexes (p=0.000), while only among females was L h significantly 
greater for wild frogs (p=0.000). L h /M b was used as an index for origin 
discrimination to eliminate the influence of M b due to variation of hus¬ 
bandry conditions among farms. Mean L h /M b for farmed frogs was sig¬ 
nificantly lower than for wild frogs (p=0.000) in the 2-year old age class. 
Discrimination correctly classified 84.4% of farmed and 96.3% of wild 
male frogs. Among females, 92.9% of farmed frogs and 90.1% wild frogs 
were correctly classified. The overall correctness of classification was 
92.1% and 90.8% for males and females, respectively. However, L h /M b 
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revealed variation with age, resulting in reduced discriminative power for 
frogs >3 years old. We introduced a coefficient C e to adjust the L h /M b of 
frogs >3 years to the level equivalent to 2-year frogs. The adjustment 
achieved 89.5% for overall correctness of origin for wild males and 
92.4% for wild females >3 years old. These results show that L h /M b is an 
effective index to discriminate wild from farmed Dybowski’s frog. Since 
the physical demands of jumping are common among anurans, this index 
is also potentially applicable to other anuran species. 

Key words: Dybowski’s frog; morphometric discrimination; wild, 
farmed;, Rana dybowskii 

Introduction 

Humans have always eaten wild animals. However, wildlife has 
been depleted worldwide and is no longer adequately abundant 
to serve as a human food resource in most of the world. To com¬ 
pensate for this shortage, wildlife farming has been adopted for 
some high-value species. Many wildlife species have been do¬ 
mesticated and farmed successfully, including fur-bearing mam¬ 
mals, snakes, frogs, crocodiles, ostriches, deer, giant rat, snail, 
geese, and ducks. Many farmed wildlife species are endangered 
in the wild and protected by law or regulation. However, their 
wild populations are still occasionally harvested or poached and 
these wild animals enter food markets alongside their farmed 
conspecifics. This raised the question of how to distinguish the 
wild animals from those reared on farms. 

One approach to discrimination of wild from farmed animals 
is based on differences between the wild and farm environments 
and the impacts of these differences on animal physiology. It is 
widely accepted that morphology and physiological processes 
evolve in response to environmental stimuli to enhance individ¬ 
ual survival and reproductive success. In contrast to nature, the 
captive environment is confined, requires lower activity levels, is 
limited in terms of diet, requires less energy expenditure to pro¬ 
cure food, and has a lower predation risk (McPhee 2004). Cap¬ 
tive animals respond to this less demanding environment through 
physiological and morphological evolution that results over time 
in phenotypic divergence from their wild conspecifics (McPhee 
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2003a; McPhee 2003b). These evolutionary changes can point 
researcher and law enforcement personnel toward better ap¬ 
proaches to segregating wild from farmed animals. 

Some researchers developed analytic methods based on stable 
isotope deposition in hard tissues such as fur, claws, and teeth 
(Hammershoj et al 2004). Carbon ratios ( 13 C/ 12 C) have been used 
to distinguish farmed from wild animals based on differences in 
their food sources (Chisholm et al 1982; Kelly 2000; Smith et al 
1996; Knoff et al 2002). A similar approach is based on the dif¬ 
ferences in fatty acid profiles between wild and captive animals 
due to differences in sources of dietary fat (Rouvinen et al 1989). 
However, both stable isotope content and fatty acid profiles may 
change in response to changes in diet (Ostwald 1962; Knoff et al 
2002), thus are only applicable for special situations. Bone struc¬ 
ture and density are associated with diet and activity. Animals in 
the wild tend to be more active than those in the farm; therefore, 
femur bone density is a stable index useful for discriminating 
wild and farmed animals (Yang et al. 2011). However, these ap¬ 
proaches require bones to be isolated from the test animals and 
are often impracticable for use on live animals. 

The provision of an adequate, readily obtainable diet together 
with protection from potential predators and provision of con¬ 
fined living spaces would cause captive animals to spend 
significantly less energy than wild animals on feeding, foraging, 
traveling and vigilance. As a result, the body fat mass of captive 
animals should be higher, while activity and motion intensity 
should be lower than among wild animals. 

For frogs, jump distance and frequency is determined by the 
structure and function of the hindlimbs (Peplowski et al. 1973). 
Jumping is a physical factor stimulating the development of 
hindlimbs of frogs (Olson et al. 1998.) because it positively im¬ 
pacts development of bone and muscle (Watanabe et al. 2007; 
Gary et al. 2001). As a result, hindlimbs of wild frogs can be 
considered optimally developed to provide leaping power needed 
for survival in the wild (Essner et al 2010; Leisler 1977). This 
suggests that metrics of the hindlimbs of frogs could provide a 
novel approach to discriminate wild from farmed frogs because 
frogs living in the wild could have longer hindlimbs than their 
counterparts who have lived many generations in farms. 

Dybowski’s frog (Rana dybowskii ) occurs in northeast China, 
Mongolia, Russian Far East, Korean Peninsula, and Japan (Xie et 
al. 1999). The extensive use by humans of dried frog oviduct as a 
natural invigorant causes unsustainable taking of Dybowski’s 
frog from the wild and this has resulted in declines in wild popu¬ 
lations (Liu et al. 2007). As wild populations declined, market 
prices rose, and this stimulated captive breeding trials (Li et al. 
1987; Meng et al. 1990). Dybowski’s frog farms were estab¬ 
lished after the mid 1980s throughout the species range in north¬ 
east China (Liu et al. 2007). Frogs are commonly kept in pens as 
small as 200-1000 m 2 and surrounded by the fences (Bai et al. 
2006.). They are fed a limited range of readily available inverte¬ 
brate species, such as yellow meal worm, screwworm, and 
earthworm. We hypothesized that these living conditions would 
result in changes in the morphology of frogs, specifically a re¬ 
duction in the ratio of hindlimb length to body mass, that would 
enable discrimination of farmed from wild frogs. We used Dy- 
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bowski’s frog to test this hypothesis. 

Materials and methods 

A total of 335 Dybowski’s frogs were captured in the Jiaohe 
(126°53 / E, 44°02'N, altitude 244 m), Fusong (127°33' E, 42°21' 
N, altitude 780m), and Tonghua (126° 12' E, 42°03' N, altitude 
632 m) regions of China’s Changbai Mountains in October 2006, 
including 182 females and 173 males. Age of each frog was de¬ 
termined by the tibiofibula slice procedure as described in Yang 
et al (2011). 60 frogs (32 males and 28 females) at the age of two 
years were collected from a frog farm located in the same region 
(126°29' E, 44° 10' N, altitude 174 m) where the wild frogs were 
captured. The frogs were raised in 5 mxlO m pens after meta¬ 
morphosis with density varying from 100-200 frogs per square 
meter. The frogs were harvested before winter hibernation at the 
age of two years, usually during September and October. 

Body mass (M b ) of these frogs was measured using an elec¬ 
tronic balance with 0.0lg sensitivity (Shanghai Yueping Instru¬ 
ment Co. Ltd., China). Hindlimb length (L h ) was measured from 
the hind end of the body to the tip of the fourth toe. L h measure¬ 
ments were preformed using a screw thread micrometer calliper 
with accuracy of 0.1mm. 

A relative index, L h /M b , was computed for each sex and age 
group of wild and farmed frogs. Variation of L h /M b with age and 
geographic location was examined using ANOVA multiple com¬ 
parisons (Tamhane method). Mean L h /M b values were compared 
between farmed and each of three wild populations of 2-year old 
frogs for each sex using ANOVA multiple comparisons (Tam¬ 
hane method). We further pooled the data of three wild popula¬ 
tions and compared L h /M b as a whole population to farmed popu¬ 
lation for each sex using nonparametric Mann-Whitney U tests 
(two-tailed) at significance level a=0.05. 

In order to improve the power of discrimination for frogs older 
than two years, the relationship between age and M b was esti¬ 
mated for the pooled wild populations for each sex using the 
curve estimation method of regression. The regression function 
was used to estimate the age of a given wild frog by using its 
value of M b . A coefficient C e was calculated as the mean L h /M b 
of 2 year frogs divided by mean L h /M b of 3 year, 4 year and >5 
year frogs. Then, an adjusted L h /M b of the frog > 3-year old was 
deduced through the observed L h /M b value multiplied by C e . The 
deduced L h /M b was then used for origin discrimination. Finally, 
we tested the effectiveness of the deduced L h /M b in discrimina¬ 
tion of wild frogs > 3-year old from 2-year old farmed frogs 
using discrimination analysis. All statistical analyses were per¬ 
formed using SPSS 16.0 software (SPSS Inc., Chicago, IL, USA). 

Results and discussion 

Origin discrimination of 2-year old frogs 

We collected 172 two-year old wild frogs, 52 in Tonghua includ¬ 
ing 23 males and 29 females, 59 in Jiaohe including 29 males 
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and 30 females, and 61 in Fusong including 29 males and 32 
females. The number of 3-year old wild frogs was 60 including 
30 males and 30 females. The number of 4-year old frogs was 61 
including 30 males and 31 females. Frogs > 5-year old included 
32 males and 30 females, 62 individuals in total. 

M b of male farm frogs was 16.33 g on average with a standard 
deviation (SD) 3.27. Mean M b of male wild frogs (pooled) was 
11.65g with a SD 2.52, significantly lower than farmed individu¬ 
als (p=0.000). M b of female farm frogs was 20.55±3.98 g. M b of 
female wild frogs (pooled) was 13.04±3.55 g, also significantly 
lower than farm frogs (p=0.000). L h was 51.21=1=3.49 mm in male 
farm frogs and 51.04=1=3.55 mm in male wild frogs, and this dif¬ 
ference was not significant (p=0. 791). Among females, L h was 
55.40±4.61 mm for farm frogs and 51.99±2.95 mm for wild 
frogs, a highly significant difference (p=0.000). 

Means and standard deviations of the ratio L h /M b for wild and 
farmed frogs are shown in Table 1. L h /M b exhibited geographical 
variations in both male and female wild frogs. Among males, 
L h /M b did not differ between the Jiaohe and Fusong populations 
(P=0.989). However, L h /M b for the Tonghua population was 
significantly greater than for the Jiaohe and Fusong populations 
(p=0.000). Such geographical variations were also seen among 




Lh/Mb 


Fig. 1 Frequency distribution of L h /M b i 


females. The Fusong and Tonghua female frogs had similar 
L h /M b ratios (p=0.609). However, Jiaohe females had signifi¬ 
cantly lower L h /M b ratios than did females at Fusong and 
Tonghua (p=0.000). 


Table 1 Average and standard deviation of L h /M b in wild and farmed 
frogs aged 2 years 


Origin 

Male 

Female 

Wild populations 



Jiaohe 

4.139±0.518 

3.299±0.267 

Fusong 

4.231±0.575 

4.669±0.704 

Tonghua 

5.524±1.057 

4.964±0.963 

Pooled 

4.565±0.942 

4.312±1.003 

Farmed population 

3.099±0.406 

2.746±0.377 


Comparison of farmed with wild frogs showed that for both 
sexes, mean L h /M b for farmed frogs was significantly lower in all 
three wild populations (p=0.000). A similar result was obtained 
when data of three wild populations were pooled (p=0.000). Fig. 
1 shows the frequency distribution of indexes between farmed 
and wild (pooled) frogs. 




Lh^lb 

2-year farmed and wild Dybowski’s frogs 


Discrimination of 2-year old males using L h /M b correctly clas- females, 92.9% of farmed frogs and 90.1% wild frogs were cor- 

sified 84.4% of farmed frogs and 96.3% of wild frogs. Among rectly classified. The overall correctness of classification was 
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92.1% and 90.8% for males and females, respectively. The 
Fisher’s linear discriminant function for male farmed frogs was 

Y x =4.519 LjM b - 8.264 

and for male wild frogs was 

Y 2 =6.656L h /M b -15.526 

The function for female farmed frogs was 

Y 1 =3.404LjM b -6.121 

and for female wild frogs was 

Y 2 =5.344 LjM b -11.789 

A frog could be classified as farmed if Y x > Y 2 , otherwise it 
would be classified as wild. 



Origin discrimination of >3 -year old frogs 

Variation of L h /M b with age of wild frogs is shown in Figure 2. 
In both male and female groups, L h /M b tended to decline with the 
age increment. Discrimination power between 2-year old farmed 
frogs and wild frogs older than 2 years resulted in 0 correctness 
for both male and female farmed frogs. M b of frogs at the ages of 
3 years, 4 years, and >5 years were similar (Figure 3). Therefore, 
frogs older than 2 years could be pooled to simplify discrimina¬ 
tion. Pooling yielded two age classes for wild frogs, namely, 2 
years and >3 years. We introduced an adjustment procedure: 
firstly, to assign frogs to an age class using M b data; secondly to 
create a coefficient by which L b /M b values for frogs >3 years 
could be adjusted to levels equivalent to those of the 2-year age 
group; and thirdly to discriminate the origin of frogs using 
Fisher’s linear discriminant functions. 



Fig. 2 Plots of L h /M b against age classes in wild Dybowski’s frogs 




Fig. 3 Plots of M b against age classes in wild Dybowski’s frogs 


We tested the effectiveness of the above procedure using our 
data for wild frogs. Assignment to age class can be performed 
through the discriminant functions for male wild frogs: 


Y x =l.U6M h - 7.366 

for group of 2 years, and 
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Y 2 = \.162M b -16.476 

for group of >3 years. The functions for females are: 

Y x = 0.398M,- 3.288 

for group of 2 years, and 

Y 2 =0.82 \M b -11.720 

for group of >3 years. An individual could be assigned to 2-year 
old group if Y1 > Y2, and otherwise to >3 years group. 

The equivalent adjusting coefficient C e of L h /M b was calcu¬ 
lated as the ratio of the mean L h /M b of the 2-year group to the 
mean L b /M b of the >3-year group. The coefficient, C e , we calcu¬ 
lated was 1.382 for male frogs and 1.797 for female frogs. 

Discriminant analysis showed 82.7% of 2-year old male frogs 
and 80.4% of >3 years old male frogs could be correctly assigned 
to age classes, with an overall correct percentage of 81.5%. The 
correct percentage for female assignment was 94.5% for 2-year 
old frogs, and 85.7% for >3 years old frogs, overall 90.1%. By 
adjusting with C e , the mean L b /M b of male wild frogs >3-years 
old increased to 4.566±0.701 from 3.304±0.508, while the mean 
L b /M b of female frogs increased to 4.262±0.908 from 
2.371±0.505. 

Discrimination using the adjusted L b /M b values correctly clas¬ 
sified 75.0% of 2-year old male farmed frogs and 94.6% of male 
wild >3-year frogs. Among females, 89.3% of 2-year old farmed 
frogs and 93.4% of wild >3-year frogs were correctly classified. 
The overall correctness of classification was 89.5% and 92.4% 
for males and females, respectively. 

Discussion 

Morphometric measurement does not injure the animals, thus is a 
convenient method to discriminate wild from farmed frogs. 
However, discriminating wild from farmed frogs as presented 
here is based on the measurement of morphological traits that are 
susceptible to modification by environmental conditions. Such 
traits must meet several criteria: (1) they must be in response to 
specific factors that are substantially different between captive 
and wild environments; (2) they must adequately plastic to per¬ 
mit partitioning of wild from captive specimens; and, (3) they 
must be characterized by low within-group variation so as to 
achieve high discriminative power in statistics. 

Our hypothesis was that wild frogs would have extended 
hindlimbs. However, this hypothesis was poorly supported by 
our data for hindlimb length (L b ). Dybowski’s frog demonstrated 
significant dimorphism that L b was significantly greater for wild 
than farmed female frogs (p=0.000), while the difference was not 
significant (p=0.791) among males. 

For males, M b of farmed frogs was significantly greater than 
for wild frogs (p=0.000). For females, M b was also significantly 
greater for farmed than wild frogs (p=0.000). 

Considering the fact that the load on limbs during jumping and 
landing is associated with body mass, we used M b as an objective 
reference to obtain a relative index L b /M b . Our results demon¬ 


strated L b /M b exhibited geographical variations in wild frogs. 
This demonstrates that frog growth is influenced by habitat con¬ 
ditions, and therefore absolute indexes might not be reliable indi¬ 
cators of frog origins. Despite the geographical variations found 
among the three sampled wild populations, L b /M b in 2-year old 
farmed frogs was significantly lower than that in each of the 
three wild populations at the same age for both sexes (p=0.000). 
A similar result was obtained when data for the three wild popu¬ 
lations were pooled (p=0.000). These results suggest the index 
L b /M b is a useful indicator of origin for Dybowski’s frog. Dis¬ 
crimination analysis showed an overall correctness of origin 
assignment of 2-year old frogs using L b /M b of 92.1% and 90.8% 
for males and females, respectively. We provided discriminant 
functions to perform the identification. 

However, L b /M b revealed a notable variation with age (Figure 
2), resulting in dramatically reduced discriminative power for 
frogs >3 years old. We introduced a coefficient C e to adjust the 
L b /M b values of group of >3 years to the level equivalent to 
group of 2 years. This adjustment achieved 89.5% for overall 
correctness of origin identification for wild males and 92.4% for 
wild females >3 years old. 

In conclusion, L b /M b is an effective index to discriminate wild 
from farmed Dybowski’s frog. It should be noted that most an- 
urans move by jumping and follow the same principle of 
hindlimb development. Therefore, L b /M b might also be applica¬ 
ble to similar investigations of other anurans. 
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